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Joining radio with X-rays: A revised model for SN 1993J 

C.-I. Bjornsson^ 
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ABSTRACT 

A joint analysis is done of the radio and X-ray observations of SN 1993J. It 
is argued that neither synchrotron cooling behind the forward shock nor thermal 
cooling behind the reverse shock is supported by observations. In order for adi¬ 
abatic models to be consistent, a reinterpretation of the radius of the spatially 
resolved VLBI-source is needed during the hrst few hundred days. Instead of re¬ 
flecting the position of the forward shock, it is then associated with the expansion 
of the Rayleigh-Taylor unstable region emanating from the contact discontinu¬ 
ity. Although observations imply a constant ratio between the energy densities in 
magnetic helds and relativistic electrons, they do not appear to scale individually 
with the thermal energy density behind the forward shock; rather, in adiabatic 
models, the evolution of the magnetic held strength is best understood as scaling 
inversely with the supernova radius. 

Subject headings: radiation mechanisms: non-thermal — stars: mass-loss — 
supernovae: general — supernovae: individual (SN 1993J) 


Introduction 


SN 1993J is one of the most well-observed supernovae, detected in the radio and all 
the way to soft 7 -rays. In addition, VLBI-observations spatially resolved the radio source 
early on. The interaction of the supernova ejecta and the circumstellar medium, due to the 
wind from progenitor star, creates s hock waves , which are usually taken to be described by 
the self-similar solutions derived bv IChevalierl (jl982aj). The standard model ascribes both 
the radi o and X-ray emi ssion to come from the region in between the forward and reverse 
shocks flChevalieii Il982bl) . It is therefore expected that an analysis of the radio and X-ray 
observations should give complementary information regarding the physical conditions in 
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the shocked region. Although there have been a number of papers addressing the physical 


Suzu 

d et ah 1993: Suzuki & h 

omoto 1995; 

Fransson & Biornsson 1998; Perez-Torres et ah 

2001 

Mioduszewski et ah 

2001 

, only a few 

have attempted the synergy that a joint analysis 

could give (^e.g., 

Fransson et ah 

19961. 


Many of the observed properties of SN 1993J stand out in comparison with other su¬ 
pernovae. Although the observa tions strongly indicate that the initial X-ray emission was 
domi nated by free-free emission (jZimmermann et ahlh994j: iTanaka et al.lll994 iLeising et ah 
199411 . it proved hard to model the observations with a single component evolving with time; 
for example, th e soft X-ray declined quite slowly while the harder X-ray showed a more stan¬ 
dard behavior fjChandra et al.l 1200911 . Intriguing as this may be, the most distinct features 
show up in the radio regime. 

The VLBI-observations revealed a conspic uous dynamical change at around 300 days 
(e.g., IBartel et al.ll2002l: iBietenholz et al.ll201l[l . when the expansion velocity changed from 
being nearly constant to a rapidly decreasing phase. The spatially resolved observations 
also showed that the brightness temperature was substantially below that expected for a 
homogeneous synchrotron source in which magnetic helds and relativistic electrons were 
in equipartition. Associated with the changing dynamical behavior was a distinct change 
in the variation of the synchrotron self-absorption frequency. During the initial, almost 
constant velocity phase, the value of the synchrotron self-absorption frequency decreased 
slowly before transiting to a more rapid decline. This later phase is in line with expectations 
from the standard model and similar to that in many other radio supernovae. If the evolution 
of the self-absorption frequency is related to the dynamics, spatially resolved observations 
should not be needed in order to detect such a behavior. Other well-observed, but spatially 
unresolved, radio supernovae have not given any clear indications of a similar behavior, which 
makes SN I993J quite unique in this respect. 

There has been a substantial amount of observations p ublished after most o f the models 
for SN 1993J were propos ed; in particula r , in th e radio by IWeiler et al.l (120071) and a com¬ 
pilation of X-ray data bv IChandra et al.l (120091) . The hrst aim of the present paper is to 
critically discuss some aspects of these early models to see how they fare in a comparison 
with the later observations. Although the starting point is the radio properties of SN I993J, 
it is shown in Section|2] that a joint discussion of the radio and X-ray observations provides 
important insights to a physical understanding of the shocked region. It is concluded that 
the earlier models need modifications and/or additions in order to account simultaneously 
for both the radio and X-ray emission. Section[3] presents a toy model, which remedies some 
of these shortcomings. It is argued in SectionH] that the Rayleigh-Taylor instability at the 
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contact discontinuity can give a physical content to the toy model and, hence, make it a vi¬ 
able model for the radio as well as the X-ray emission. It is also shown how such a model can 
provide a framework in which several of the observed properties of SN 1993J, not discussed 
in detail in this paper, can be understood. A summary of the main conclusions follows in 
Section[5l 


2. Modeling the radio observation of SN 1993J 


A homogeneous, spherically symmetric, expanding synchrotron source is often suffi¬ 
cient to adequately model the radio observations of supernovae. The detailed observations 
of SN 1993J, including angularly resolve d VLBI measurements fe.g.. iMarcaide et al.l 119971: 
Bartel et al.ll2002l: iBietenholz et al.ll2003[) . made it possible to test several of the simplifying 
assumptions included in the standard model. The radio emission is thought to originate from 
behind the forward shock. For most radio supernovae only two independent observables can 
be obtained (in addition to the spectral index) from observations, namely, the synchrotron 
self-absorption frequency (z/abs) and the corresponding spectral flux (F(z/abs))- In order to 
deduce the magnetic held strength and source radius some assumption needs to be made 
regarding the energy density of the relativistic electrons with respect to that in the magnetic 
helds. 

The angularly resolved observations provided a thir d observable for SN 199 3J. Although 
the assumed sphericity was conhrmed by observations flBietenholz et al.l 1200311 . the derived 
brightness temperature was lower than expected for equipartition between magnetic helds 
and relativistic electrons. In a homogeneous source, this necessitates the energy density 
of the magnetic helds to be muc h larger than that i n rela tivistic electrons. As shown by 
Fransson &: Biornsson fjl998 ') and Perez-Torres et ah ( 2001 1 the magnetic held implied is so 
strong that synchrotron cooling ahects the value and temporal evolution of z/abs- This then 
also accounts for the initially unusually slow temporal decrease of z/abs- The reason is that 
the importance of cooling decreases with time and, hence, the column density of relativistic 
electrons decreases less rapidly than expected in the standard model. 

Although the standard model with a large deviation from equipartition between mag¬ 
netic helds and relativistic electrons gives a good ht to the observations, there are a few 
implications which may suggest the need to revise some of its basic tenets. Since thes e 
were not elaborated on in either iFransson &: BiornssonI fjl998[) or iPerez- Torres et al.l (120011) , 
they are discussed in turn be low. When a quantitative comparison is needed, the results in 
Fransson &: BiornssonI (119981) will be used. 
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For a steady wind from the progenitor star, the thermal energy density behind the 
forward shock varies with time {t) as t~'^. It is often assnmed that the energy density in the 
magnetic helds scales wit h the thermal energy d ensity so that the magnetic held strength 
B oc t~^. In SN 1993J, iFransson fc BiornssonI fll998l) dednced that the magnetic energy 
density normalized to the thermal energy density is ee ~ 0.14, while the corresponding energy 
density of relativistic electrons (ce) was several hundred times smaller. An alternative scaling 
of the magnetic held corresponds to the situation when its energy density varies inversely 
with the radiating surface, i.e., for a spherically symmetric source B oc R~^, where R is its 
radius. Such a scaling implies a non-constant value for cb. 

The mechanism responsible for the amplihcation of the magnetic held behind the shock 
is not well understood. A correct description of its scaling relation based on observations is 
likely to constrain possible scenarios; for example, in order to distinguish between the t~^ 
and R~^ scaling relations, a rapidly (d e creasin g velocity of the forward shock is helpful. In 
the self-similar solutions of IChevalierl (1l982ajl . this corresponds to a low value of n (i? oc 
where n is the power law index of the density structure of the supernova ejecta. 
The angularly resolved observations of SN 1993J shows that after an initial phase of almost 
constant expansion velocity, the outer radius of the radio source makes a distinct transition 
to a strongly decelerating phase (n ~ 6 — 7) after a few hundred days. Hence, the late time 
evolution of the radio emission of SN 1993J ohers a possibility t o throw light on the scaling o f 
the magnetic field strength. Unfortunately, in the modeling bv IFransson fc BjornssonI (119981 ) 
the error bars for the H-values deduced from a procedure increase substantially 

in the decelerating late phase as compared to the early, nearly constant velocity phase. 
Although, formally, a R~^ scaling is preferred over a t~^ scaling, the error bars are large 
enough that the latter cannot be excluded. The reason for this increase of the error bars can 
be understood as follows. 


When synchrotron cooling is important, the optically thin spectral flux can be written 


vF{v) oc R^VshUeil), 


( 1 ) 


Where Ugh is the velocity of the forward shock in the standard model and 74(7) is the energy 
density of relativistic electrons with Lorentz factor 7. The value of 7 to be used in equation 
(fTj) is that corresponding to emission at v (i.e., 7 oc {v/BY^). For a distribution of electron 
energies given by N^f) oc 7“^ for 7 > 7min and p > 2, 74(7) ^ 74 ( 7 / 7 rnin)^~ ^. where 74 is 
the total energy density in relativistic electrons. IFransson fc BiornssonI fll998l) used p = 2.1 
so that the variation of 74(7) with 7 is small; in order to simplify the notation, 74(7) ~ 74 
will be used below. Since R ~ Ught, 


uF{u) (X vl^t^Ue. 


( 2 ) 














5 


This shows that in the early, nearly cons tant velocity phase, the observed flat light curves 
imply f/g oc t~‘^. In the decelerating phase, IWeiler et ahl (120071) And that all the well observed 
optically thin light curves have Fi^u) oc up to approximately 3100 days (when a rapid 
achromatic decline of all the light curves sets in). For n ^ 6 —7, this again implies f/g oc t~‘^. 
Although it is not possible to distinguish between a R~^ and t~^ scaling relation during 
the nearly constant velocity phase, the decelerating phase clearly shows that in a cooling 
scenario the energy density of relativistic electrons scales with the thermal energy behind 
the shock. Furthermore, since the variation of the energy density of relativistic electrons is 
a direct reflection of the observ ed light curves, the small error bars for the former deduced 


m 


Fransson fc BiornssonI (119981) is due to the small scatter in the light curves. 


The light curves are independent of the magnetic held strength as long as synchrotron 
cooling is important. Instead, the variation of B can be deduced from the changing values 
of ^abs 

z/fbs OC t/gf/BArgool, (3) 

where p = 2 has been used. Here, Ub = B'^/Sn is the energy density in the magnetic held 
and Argooi is the distance electrons move behind the shock before they lose most of their 
energy. The cooling time is oc B~‘^'j~^, which leads to 


■^abs 


OC 


Ue^sh 

T^abs 


Since 7abs is the Lorentz factor corresponding to i^abs, 

Z/abs OC (UeVsh)'^^'^B^^'^. 


( 4 ) 


( 5 ) 


During the nearly constant velocity phase, iFransson fc BiornssonI (119981) And that z^abs oc 
t-o-68 gives a good representation of the time dependence of the self-absorption frequency. 
With Ue oc t~‘^, as obtained independently from the light curves, B oc t~^ is implied. Again, 
the small error bars for the values of B reflect the small scatter in the deduced values of z^abs- 


For the deceleration phase, IFransson &: BiornssonI (119981) And z^abs oc However, 

in contrast to the earlier phase, the error bars on the deduced H-values are substantially 
larger. This is not due to an increased observed scatter in the values of z^abs, since it is 
similar in the two phases (the dispersion in the exponents is ±0.03). Since synchrotron 
cooling affects the value of z^abs also in the beginning of the deceleration phase, the change 
in its time dependence is due mainly to the changing velocity. It is seen from equation ([S]) 
that z^abs oc for B oc t~^ and z^abs oc vH'^ for B oc R~^. The changing behavior of z^abs 
is consistet with B oc (e.g., Ush oc for n = 6), while B oc R~^ gives a too small of 

a change. This need for a rapidly decreasing H-value is in tension with the optically thin 
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spectra and light curves, which indicate that cooling is more important than a B oc t ^ 
scaling would suggest; hence, the increased error bars for B. 

The need f or a s caling relation of B slower than t~^ is reinforced by the observations 
of lWeiler et al.l (120071) . which show a steady power law d ecline of all the optically thin light 
curves from >100 days until Ri 3100 days. In the model of iFransson fc BiornssonI fjl998l) . this 
means that if cooling was important at the beginning of the decelerating phase, it must be 
so also at 3100 days. The longest wavelength with a well observed light curve is A = 20 cm. 
This, then, gives a lower limit to the magnetic held strength, which corresponds, roughly, to 
a R~^ scaling from the nearly constant velocity phase and, hence, cb ^ 0.78. In fact, a value 
larger than this lower limit is suggested by the variation of F(i/abs)- Together with Ue oc 
equation (|2]) yields 


B{l^ahs) OC 


dL 

^abs 


( 6 ) 


which shows that in the model of iFransson fc BiornssonI fll998l) . F(z/abs) is expected to stay 
roughly con stant fi.e., indepen dent of r'abs)- This is indeed the case for the light curves 


observed bv lWeiler et al.l (120071) during the deceleration phase. It should be noticed that the 
roughly constant value of F"(r'abs) is model dependent, since it results from the combination of 
a specihc velocity law and the assumption of cooling. Although the light curve for A = 90 cm 
is not of the same quality as those at shorter wavelengths, it is consistent with a constant 
value of A(i/abs)- The maximum occurs at f ~ 2000 days, which results in an unrealistically 
large value of B (i.e., cb > 1). Such a large value of cb could be avoided by invoking a 
larger mass-loss rate from the progenitor star. This would increase the mass and energy 
requirements for the ejecta (see Section[3]). Furthermore, in order to accommodate the free- 
free absorption, a heating mechanisr n of the circurn s tehar medium has to be found that is 
more efficient than that discussed in iFran.sson et abl (119961) . 


2.1. Models with negligible synchrotron cooling 


A revised model for SN 1993J will be discussed in Section[3l for which synchrotron 
cooling is not important. The light curves in the optically thin phase are then given by 
F{u) ^ T(z/abs)(i^abs/i^)*'^~^^^^- Assumiug no cooling. iFransson fc BiornssonI (119981) deduced 
p = 2.7. Together with A(z/abs) ~ constant and i^abs cx: this shows t hat adiabatic 

model s can account for the optically thin light curv es (F(z/) f ound bv IWeiler et al. 


(120071) in the deceleration phase. Furthermore, since iBietenholz et al.l (120111 ) found 

for the outer radius of the radio source (i.e., n ^ 7), I'abs^ R~^, which is consistent with 

B oc R~^ but not B oc 
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For a standard synchrotron model, -F(i/abs) oc The rapid decline of the 

velocity in SN 1993J is particnlarly nsefnl to estimate how much the scaling of can deviate 
from that of oc for example, 74 oc t~‘^ implies -F(i/abs) oc oc 

where i^abs oc R~^ has been used. With n ~ 7 and p = 2.7, F(z/abs) SS ^'abs'^- This implies that 
between A = 3.6 cm and A = 20 c m, Ffz/ahs) should cl echne by a factor 1.9, which is too fast 
to be consistent with observation flWeiler et al.l 120071) . Hence, this suggests that also for Ue 
the scaling is closer to R~^ than t~^. A more detailed discussion of the scaling relations for 
Ue and Ub in radio supernovae will be given in a forthcoming paper. 


2.2. Constraints on the properties of the supernova ejecta and the wind from 

the progenitor star 


With a constant mass-loss rate (M) from the progenitor star, the supernova shock has 
swept up a circumstellar mass (Ms^) at time t corresponding to 


= 


n 


2 "t^sh 

- Mt—. 

n-3 Uw' 


(7) 


where Veh/v^ is the ratio of the velocities of the forward shock and the wind of the progenitor 
star, respectively. In the thin shell approximation, the ejecta mass behind the reverse shock 
is Mgj = {{n — 4)/2)Msu, which can be written 


,9 — 4)(n — 2) ,• Ugh 4 

Mej,0 = 10-2^— 


2(n - 3) 


^^w,l 


( 8 ) 


where Mgj^Q is the ejecta mass behind the reverse shock in solar mass units, M _5 is the 
mass-loss rate of the progenitor star in units of 10 “® solar masses per year, ty^ the time since 
the supernova explosion measured in years, nsh ,4 = Ush/lO^kms“^, and n*,! = nw/ 10 kms“^. 


The best-£t parameters obtained in iFransson &: BiornssonI fll998l) are M _5 = 5, t = 


100days, and nsh ,4 = 2.2. From equation (jS]), this yields = 1.5 x 10 “^(? 7 , — 4)(n — 

2 )/(? 7 , — 3). As deduced from the spatially resolved VLBI-observations, the velocity of the 
outer extent of t he radio emission mov es with roughly constant velocity, which implies a 
large value for n. IFransson et al.l fjl996[) argued that the X-ray observations from SN 1993J 


are best accounted for by n 25, which then corresponds to Mej,© = 0.33. For n 20 — 30, 
Ush = 2.2 X lO^kms”^ corresponds to an ejecta velocity Uej ~ 2.0 x 10^kms“^. The value 
of Mej ^0 deduced from equation fISll i s then more than an order of magnitude larger than 


that obtained bv lWooslev et al.l fll994l) for their best fit model (13B; O.OI 6 M 0 moving faster 
than 2.0 x 10^kms“^) to the optical observations of SN 1993J. This discrepancy between 
the deduced values of the ejecta mass is likely to be even larger, since the roughly constant 

















veloci ty phase appears to last a factor of a few longer than assumed in iFransson fc Biornsson 
fll998l) flBartel et al.ll2002l) . Other models of the supernova explosion do not seem to be able 
to bri dge this gap in ejecta masses; on the contrary, models discussed in ISuzuki fc Nomoto 
f l995h have ei e cta m asses at high velocities substantially smaller than those obtained by 
Woosley et ah (1994). 


The implications of large ejecta masses at high velocities can also be illustrated by 
considering the corresponding energy. The energy contained in the ejecta at velocities larger 


than Vej is 


-E'ej {Uei ] 


n — 3 , , o 

^--MejUp:. 

2(n -5) ^ "J 


(9) 


With n = 25 and t = 300 days flBartel et ahl l2002l) . equation ([9]) gives = 2.0 x 


lO^km/s) ~ 4.4 X 10^^ erg/s. Such large values of the energy at high velocities lie outside 
the range thought possible in standar d core collapse superno vae and, qualitatively, are more 
similar to engine driven supernovae (jSoderberg et al.ll2010l) . In order to obtain physically 
more realistic values for Mej and one needs to lower the values of M and/or n. 


Fransson et ahl (119961) considered two qualitatively different models for the X-ray emis¬ 


sion from SN 1993J. One model assumed cooling to be important behind the reverse shock. 
This results in absorption of the X-ray emission by the cold gas and, hence, in itially no X-ray 


emissi on is expected to be observed from the reverse shock. As emphasized in iFransson et ah 


fll996l) . the most characteristic aspect of this radiative scenario is the evolution of the X-ray 
flux. After about a hundred days, the cold gas becomes transparent and the reverse shock 
starts to contribute to the observed X-ray emission. The declining X-ray light curves should 
then reverse and start to rise (ROSAT-band) or flatten out (ASCA-band). In the other 
model, cooling is assumed to be unimportant behind the reverse shock. Therefore, in this 
model, the observed flux should come from the forward as well as the reverse shock also 
during the initial phase of shock expansion. The two models imply quite different values for 
M _5 and n. In the radiative case M_r ~ 3 — 5 and n 25 — 30, while in the adiabatic 


case M _5 r:! 1 — 3 and n 6 — 8. iFransson et ahl (Il996h argued that o bserv ations favored 


the radiative model. However, as has been discussed bv I Chandra et ahl (120091) . there are no 
indications from late time observations of a break in the X-ray light curyes neither in the 
ROSAT-band nor the ASCA-band at around a few hundred days. Hence, this supports an 
adiabatic model, which has been advocated bv ISuzuki et ahl (119931 1. and is in line with the 
conclusion drawn above. 

The density of the circumstellar medium is normally taken to be a power-law, oc 
for a steady wind, s = 2. The complexities of both the radio and X-ray emission from SN 
1993J haye sometimes been attributed to a non-steady wind, i.e., s ^ 2. ISuzuki fc Nomoto 
(I 1995 I) modeled the X-rays with s <2 during an initial phase, which transits at around day 
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40 to a phase with s > 2. This causes the velocity of the forward shock to decrease faster 
initially, while at later times slower, than for s = 2. Qualitatively, this expected flattening 
of the variation of the shock velocity wi th time is opposite the behavior observed by the 
spatially resolved VLBI-measurements of iBartel et ahl (120021) . 


It is generally agreed that synchrotron self-absorption as well as free-free absorption 


Fransson & Biornsson 

1998; 

Chevalier 

1998; 

Perez-Torres et al. 

2001; 

Mioduszewski et al. 


200ll) . Even so, this does not give a good £t for s = 2, since, for example, the resulting light 
curves are too narrow. As discussed in Section^ this can, at least in the initial phase, be 
remedied by assuming a magnetic held strong enoug h to cause synchrotron cooling. Another 
possibility suggested bv iMioduszewski et al.l (120011) is that s < 2. Since the time variation 
of the free-free absorption (rg) is deduced from observations, this has implications f or the 


temperature structure of the circumstellar medium (Tcm)- iFransson fc BiornssonI (119981) 
showed that s = 2 together with a declining value of Tcm ahead of the forward shock 
resulted in free-free absorption consistent with observations. With r g oc where 

Up PC P w is the density of thermal electrons, the value of s deduced by 
(1200 ll) instead implies Tcm ~ constant. 


Mioduszewski et ah 


In IFransson fc BiornssonI (119981) . the variation of Tcm was taken from self- consistent cal- _ 

culati ons of the Compton heating by the X-ray emission from the shock, while in IMioduszewski et al 
(120011) it was treated as a free parameter. The much lower densities in the latter model irn - 
plies negligible Compton heating. However, as discussed in iBiornsson fc LundqvistI (120141) . 
the free-free absorption itself causes heating of the circumstellar medium. This then sets 
a lower value for Tcm- During the initial phase of radio emission from SN 1993J, when 
free-free absorption is most promi nent, this lower limit is a f actor of a few larger than the 
value Tcm = 7 x 10^ K deduced by Mioduszewski et al. ( 2001 ). 


As will be discussed further in Section|H a distinguishing feature of the radio emission 
from SN 1993J is the low brightness temperature. This is another observational constraint 
that is not directly encompassed by models with s < 2, since the value of s should not affect 
the brightness temperature. Together, this shows that there is no compelling evidence for 
s 7 ^ 2 neither from X-ray nor radio observations. 


It was argued in IBiornsson fc LundqvistI (120141 ) that a self-consistent calculation of the 
temperature of the circumstellar medium, assuming heating only by the observed free-fre e 
absorption, gives a fit at least as good as the one obtained in IFransson fc BiornssonI (119981) . 
The main free parameter is the density of the circumstellar medium, i.e., the mass-loss rate 
of the progenitor star. The deduced value lies in the range M _5 = 0.8 — 1.0 for Uw,i = 1. The 
low density of the circumstellar medium implied by such a mass-loss rate makes additional 
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heating mechanisms unlikely. The VLBI-observations show the velocity, interpreted as that 
of the forward shock, after about a year, to vary as n ^ 7 . An important assumption 
underlying the claim of a good fit in iBiornsson fc Taindqvistl fj2014j) was that this velocity- 
law also applies to the initial phase, i.e., during the hrst few hundred days of shock expansion. 
If this were the case, the observed VLBI-velocity during this phase cannot correspond to that 
of the forward shock. The implications of such a model are discussed in the next section. 
It may be noticed that M _5 = 0.8 — 1.0 and n Rr! 7 give values of the ejecta mass at high 
velocities consistent with that of model 13B in 


Woos. 


values are also close to those deduced by ISuzuki et al 
few weeks of X-ray observations. 


ev et al.l (119941 ). Incidentally, these 


( 1993 !) from an analysis of the hrst 


3. A revised model for SN 1993J 


In previous modeling of the radio and X-ray emission from SN 1993J, it has tacitly 
been assumed that the velocity deduced from the expanding VLBI-source corresponds to 
that of the forward shock. In this section, a revised model is discussed in which this is not 
the case. It takes its starting point in the VLBI-observations, which show the outer radius 
of the radio source to increase, initially, almost linearly with tim e with a later transition 
to a slower expansion flBartel et al.l l2002l: iBietenholz et al.N201lh . Extrapolating the later 
expansion back to the earliest VLBI-observations, one hnds that the difference between 
the extrapolated and measured radii is approximately equal to that expected be tween the 
forward shock and contact discontinuity/reverse shock in the self-similar models of [Chevalier 
(jl982al ). A simple toy model to depict such a situation is shown in Figure[T] Since the 
distance between the contact discontinuity and the reverse shock is much smaller than the 
distance between the forward shock and the contact discontinuity, no distinction is made 
between contact discontinuity and reverse shock in the discussion below. The radii of the 
contact discontinuity (i?c) and forward shock (i?Fs), respectively, are assumed to vary with 
time according to the self-similar model. The radius R is the outer radius of the synchrotron 
source; i.e., the synchrotron emission is produced between R and Rq. The observed transition 
of the variation of R at around 300 days is assumed to correspond to R i?FS (or some 
constant fraction thereof). At later times, the expansion of the VLBI-source is determined 
by the forward shock. 


Such a model can provide alternative explanations to both the radio and X-ray obser¬ 
vations. As discussed in Section[2l the initial decrease of the self-absorption frequency (z/abs) 
with time was slower than expected for an adiabatic source. With synchrotron cooling, the 
thickness of the synchrotron emitting shell increases faster than for the adiabatic case and 
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observations could be accounted for by standard assumptions regarding the variation of the 
magnetic held strength and density of relativistic electrons. Hence, with an appropriate 
choice of R{t), the same result can be obtained for an adiabatic source model. 


The roughly constant velocity of the VLBI-source during the hrst few h undred days im¬ 


plies a large value for n in the standard, self-similar model. As discussed in iFransson et ah 


(119961) . the reverse shock is then cooling and the observed X-ray emission comes from behind 
the forward shock. In the proposed revised model, the value of n is not constrained during 
the initial phase and, for example, it may be the same as that appropriate for the later phase. 
Hence, the initial VLBI-measurements cannot be used as arguments against adiabatic mod¬ 
els. Rather, arguments in favor of adiabatic r nodels can found fr om two independent aspects 
of the observations, name ly, X-ray emission fISuzuki et al.l 119931) and free-free absorption of 
the synchrotron emission ( Biornsson fc LundqvistI 2014 ). both of which give similar values of 
M _5 and n. It may be noticed that in order for adiabatic X-ray models to be applicable, it 
is necessary for the synchrotron emitting region to grow outwards starting from the contact 
discontinuity rather than inwards starting from the forward shock. 


The rate of increase of the outer VLBI-radius changes at a few hundred days. iBartel et al 


fl2002l) show that a good £t can be obtained with two power-laws, which cross at approxi¬ 
mately 300 days. If this were due to increased deceleration of the forward shock, one would 
expect the transition to take place over at least a dynamical timescale. In such a situation, 
the measured radii during the transition phase should lie below the extrapolations of both 
power-laws. The observed error bars are small enough to argue that this may not be the 
case. If so, this would indicate a rather abrupt transition more akin to a discontinuity in 
the velocity. Such a discontinuity is straightforward to account for in the revised model but 
harder to incorporate into a model where the increasing deceleration of the forward shock is 
due to external changes in either the circumstellar matter or the supernova ejecta. 


In principle, the VLBI-observations could be used to deduce R{t). Hence, for a given 
value of n, and assuming a homogeneous distribution of relativistic electrons and magnetic 
fields in the region R — Rq, the synchrotron self-absorption could be calculated and compared 
to the observed r'abs(^)- However, there are several effects that limit the usefulness of such a 
comparison. As will be discussed in the next section, at least the magnetic held is likely to 
be highly inhomogeneous in the synchrotron emitting region. Furthermore, during the initial 
phase, the outer and inner radii of the VLBI-source could not be determined independently. 
Instead a constant ratio was assumed, which was taken to be the same as that observed during 
the later phase of expansion. Such an assumption is inconsistent with either synchrotron 
cooling models or the proposed revised model. Although this is not expected to change the 
deduced values of the outer radius substantially, it can have a non-negligible effect on z/abs(t). 
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In order to illustrate a few of the implications of the revised model, let n = 7. For 
a given observed radius, the instantaneous velocity depends on n as v oc (n — 3)/(n — 2). 
Hence, at an assumed transition on day 300, the velocity of the forward shock is smaller 
than in the nearl y constant velocity scenario by a factor 0.8, i.e., = 2.2 x 0.8. From 

Chevalierl fjl982a|l . one also finds that fej = 0.93Ugh- The variation of the ejecta velocity 
at the reverse shock is then given by Uej = 1-6 x 10^(f/300days)“°'^kms“^. The first radio 
and X-ray observations took place at t ~ 10 days, which corresponds to an ejecta velocity 
3.2 X 10^kms“^. This is close to the maximum velocity of model 13B in IWooslev et ah 




fll994j) and shows that this model is consistent with a self-similar solution also during the 
early observations. 


Weiler et al.l (120071) have emphasized that an achromatic break occurs in all the radio 
as well as X-ray light curves at t ~ 3100days, which corresponds to ~ 1.0 x 10^kms“^ 
This value matches closely the ejecta velocity in model 13B, where the density distribution 
makes a sharp break and enters a region where it stays almost constant with decreasing ra- 
dius/velocity. This region represents the inner part of the hydrogen rich envelope. A possible 
cause for the achromatic breaks in the light curves is, therefore, that the energy/momentum 
from the ejecta is no longer large enough to maintain a self-similar shock structure; i.e., 
the reverse shock would weaken or even disappear. Not only would this reduce the X-ray 
emission but also, as is discussed in the next section, cause a decli ne in the radio emission . 
Support for such a scenario comes from the VLBI-observations of iBietenholz et al.l (120111) . 
After a few years, the radio source had expanded enough for the observations to allow in¬ 
dividual measurements of the outer and inner radii of the synchrotron emitting shell. The 
outer as well as the inner radius follow a n ~ 7 evolution until ~ 3000 days. At this time, 
the expansion of the inner radius slows down considerably, while the outer radius is much 
less affected. In the revised model, the base of the radio emitting region is the contact dis¬ 
continuity. A weakening of the reverse shock is expected to affect the velocity of the contact 
discontinuity and, hence, the velocity of the inner radius of the radio source. 


4. Discussion 


The toy model shown in Figured] envisions a synchrotron emission region that starts 
to grow outwards from the contact discontinuity to fill all, or a constant part, of the re¬ 
gion behind the forward shock. Such a scenario accords well with the physical situation 
expected to emerge from the Rayleigh-Taylor instability at the contact discontinuity. Den¬ 
sity contrasts by more tha n a factor 10 result as shocked ejecta and circumstellar matter 
mix (IChevalier et al.lll992l) . Likewise, the strength of the magnetic held amplihed by the 
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turbulence varies considerably, with the highest values in regions bounding the high den - 
sity hngers of ejecta gas protruding from the contact discontinuity fl.Tnn fc Normanlll996al) . 
Roughly half of the volume between the forward and reverse shocks is estimated to be af¬ 
fected by the instability in its saturated phase. This fraction is likely to increase with the 
compression ratio across t he forward shock; for example, due to a signihcan t amount of 
relativistic electrons/ions (IBlondin fc Ellisonll200ll: iDuffell fc MacFadvenll2014f) . The phys¬ 
ical underpin ning of the toy model in S ection|3] suggested here is then similar to the one 
discussed by IChevalier et al.l (1199211 and lJun fc NormanI fjl996a|l for the radio emission in 
SN 1572 (Tycho) not directly associated with the forward shock. The association of the 
synchrotron source with the turbulent region implies two things. Firstly, the amplified mag¬ 


netic fi eld should dominate the magnetic held behind the forward shock fe.g.. lJun fc Norman 


1996bl) . Secondly, in order for the synchrotron emission region to be dehned by the magnetic 
held distribution, the relativistic electrons should hll this region more or less uniformly; for 
example, through hrst order Fermi-acceleration at the forward shock. 

The ehects of the amplihcation of the magnetic held by the R ayleigh-Taylor ins t abilit y 
on the spatially resolved VLBI-source have also been considered in iBietenholz et al.l (1200311 . 
The focus was on the evolution of SN 1993 J after a few years, when the outer and inner radii 
could be resolved individually and the source was in the strongly decelerating phase. They 
concluded that the outer radius is rather close to the forward shock but hnd no indications 
for the shock-front itself to be ahected by turbulence. In contrast to the revised model 
discussed in this paper, they suggested that the two diherent velocity regimes may indicate 
that the shock structure was not self-similar. 

The radio observations of SN 1993J clearly showed the brightness temperature to be 
below that expected for a homogeneous synchrotron source in which the magnetic field is in 
equipartition with the relativistic electrons. The slower than expected decrease of the syn- 
ch rotron self-absorption f r equen cy increased the brightness temperature with time. However, 


m 


Fransson fc BiornssonI (119981) the brightness temperature is also lowered by the dominance 


of the magnetic energy density over that in relativistic electrons, which provides a time in¬ 
dependent factor. In the revised model, the initial increase of the brightness temperature 
and the associated slow decrease of the synchrotron self-absorption frequency are due to the 
relative increase of the thickness of the synchrotron emitting shell, while the constant factor 
is attributed to the inhomogeneous str ucture of the mag netic field, which results in an effec¬ 
tive covering factor smaller than unity (jBiornssonll2013l ) . A more quantitative description of 
the variations of brightness temperature and synchrotron self-absorption frequency is harder 
to do, since this is sensitive to the details of the inhomogeneous source structure. 


As discussed in Section[3], the break in the X-ray light curves at f ~ 3100 days may be 
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due to the reverse shock entering a roughly flat portion of the density distribution of the 
ejecta, which makes the self-similar solutions inapplicable. With the revised model discussed 
above, the simultaneo us break in the r a dio ligh t curves can be ascribed to the same physical 
effect. As discussed in I.Tnn fc NormanI fjl996al) . a weakening of the reverse shock causes the 
density contrast at the contact discontinuity to decrease. This, in turn, weakens the driving 
of the Rayleigh-Taylor instability and, hence, leads to less amplification of the magnetic held. 
The accompanying increase of the region behind the reverse shock is then the likely reason 


for th e slowing down of the expansion of the inner boundary observed by iBietenholz et ah 
f)201ll) at the same time. 


Some of the aspects of the revised model are also co nducive to limit i ng th e possible 
interpretations of the X-ray observations. It was argued in iFransson et ahl (119961) that the 


low t emperature deduced for the X-ray emitting gas around 200 days fjZimmermann et ah 


19941 ) is h ard to accommodate within an adiabatic model. Howev er, later observations have 
indicated (jZimmermann fc Aschenbachl l2003l: ISwartz et al.l l2003l) that there is a range in 
temperatures. Although a co nsistent description has not yet be en achieved, the implied 
range is at least a factor of 10. Zimmermann fc Aschenbachl (2003) fitted the X-ray emission 
up to ~ 3000 days with a simple two temperature model, which indicates that the average 
temperature decreases with time but that the range is more or less constant. In light of the 
revised model, it is interesting to note that the evolution of both temperatures is consistent 
with equip artition be t ween ions and electrons and n 7. A similar conclusion has been 
reached bv lUno et al.l (120021) when fitting the ASCA observations during the first ~ 600 days. 


Such a range in temperatures may be caused by the turbulence driven by the Rayleigh- 
Taylor instability at the contact discontinuity. When the instability has saturated, it is likely 
to cover a substantial part of the shocked circumstellar medium. In the homogeneous, adia¬ 
batic model, the X-ray emission should come mainly from the shocked ejecta. However, the 
large density variations in the turbulent region could result in a non-negligible contribution 
of the shocked circumstellar medium to the X-ray emission even in the case of an adiabatic 
reverse shock. Although the relative contribution to the total X-ray emission from the two 
shocked regions is rather sensitive to the details of the revised model, there are a few as¬ 
pects that distinguishes it from homogeneous models with an initially nearly constant shock 
velocity (i.e., high value of n). 


The higher shock velocities early on in the revised model as compared to constant veloc¬ 
ity models increase the equipartition temperat ures behind the sho cks; for example, assuming 
equal numbers of hydrogen and helium atoms flBaron et al.lll9941) . the equipartition temper¬ 
ature behind the reverse shock is ~ 10® K at 10 days. Although, in order for the electrons to 
reach this temperature, mechanisms in addition to Coulomb collisions are needed to transfer 
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energy from the ions to the electrons, it shows that the reverse sho ck may have contributed 
to the high energy emission detected by OSSE flLeising et ahlll994[) . The difference between 
the equipartition temperature and the temperature reached by the electrons behind the re¬ 
verse shock, assuming Coulomb collisions only, is not that large. The opposite applies for 
the electrons behind the forward shock, where this difference is much larger. The turbulence 
driven by the density contrast at the contact discontinuity amplifies the magnetic field. This 
process may also excite plasma waves, which could contribute to the energy transfer between 
ions and electrons. If so, this would further enhance the relative importance of the forward 
shock to the X-ray emission. 

When the electron temperature (Te) is belo w its equipart ition value and determined by 
Coulomb collisions, it varies as Tg oc (ngTion^)^'^^ flSpitzerlllQGSl) . where Tion is the ion temper¬ 
ature. Since Tjon oc in a steady wind environment Tg oc f is expected. This scaling 
relation coincides with that for the equipartition temperature when n = 7. The deduced 
temperature variations discussed above may then not necessarily indicate that equipartition 
conditions apply. 

The observed X-ray emission limits the contribution from Comptonized flux. This leads 
to a tension between the high temperature indicated by the OSSE-observatio ns and the high 


mass- loss rates needed to account for the softer X-ray emission. As is seen in iFransson et ah 


fll996l) . this causes considerable strain when trying to hnd acceptable parameters in radiative 
models. This constraint is eased in adiabatic models; for example, M_^ = 1 and n = 7 give 
an optical depth to electron sc attering around 10 day s that is more than a factor 10 smaller 
than for the best £t model in IFransson et al.l (119961 1. This then allows for correspondingly 
higher temperatures. 


5. Conclusions 

The main conclusion of the present paper is that a revision is needed of the model(s) 
used to describe the observed behavior of SN 1993J. This is based on the following two main 
results from the analysis in Section^ 

1) Together, radio and X-ray observations show that there are a number of arguments 
against models invoking cooling in either spectral regime; for example, the implied mass and 
energy at high eject velocities are at least an order of magnitude larger than expected in 
standard core collapse scenarios. 

2) Adiabatic models can account for the radio as well as X-ray observations; for ex¬ 
ample, a mass-loss rate of the progenitor star around lO“®M0yr“^ (for a wind velocity 
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Vvi = lOkms and an ejecta density distribution corresponding to n ~ 7 are consistent 
with both radio and X-ray observations. 


The crucial new component in the revised model is the interpretation of the radius of 
the spatially resolved VLBI-source. 

3) During the first few hundred days, when the velocity of the radio supernova is roughly 
constant, the observed radius does not correspond to that of the forward shock. Instead, 
it is associated with the expansion of the Rayleigh-Taylor unstable region starting from the 
contact discontinuity. After saturation of the instability, the outer radius of the VLBI-source 
is determined by the forward shock. Hence, the velocity of the forward shock is determined 
by n ~ 7 also during the initial phase of evolution of the supernova. 

4) The radio emission results from the amplihcation of the magnetic held by the turbu¬ 
lence in the Rayleigh-Taylor unstable region. Such an explanation has previously been put 
forth for the radio emission sometimes observed to peak inside the forward shock in super¬ 
nova remnants, for example, SN 1572 (Tycho). Hence, SN 1993J would offer the possibility 
to study also the time evolution of this instability. 


5) A low value of n makes it easier to determine whether the value of a given quantity 
scales with radius {R) or time (f). In an adiabatic model for SN 1993J with n ^ 7, observa¬ 
tions suggest that the magnetic held {B) varies as R oc 1/R rather than B oc 1/f; i.e., the 
energy density of the magnetic held does not scale with the thermal energy density behind 
the forward shock. Furthermore, with B oc 1/R, observations indicate that also the energy 
density of relativistic electrons does not scale with the thermal energy density. 


6) Both the ejecta mass and velocity st ructure deduced for SN 1993J accord well with 
model 13R favored by IWooslev et al.l (1199411 from an analysis of the optical observations. 
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Log Time 


Fig. 1.— A toy model for the evolution of the synchrotron emission region in SN 1993J. i?FS 
and Rq are the radii of the forward shock and contact discontinuity, respectively, while R is 
the outer radius of the expanding emission region. Synchrotron radiation is assumed to come 
only from the dashed area between R and Rc- In the hgure, it is assumed that the emission 
region during the later phase extends all the way to the forward shock. As discussed in the 
text, this is not necessary, only that the outer radius R becomes a constant fraction of i?FS- 




